Geometrical enhancement of the proximity effect in quantum wires 
with extended superconducting tunnel contacts 
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We study Andreev reflection in a ballistic one-dimensional channel coupled in parallel to a su- 
perconductor via a tunnel barrier of finite length L. The dependence of the low-energy Andreev 
reflection probability Ra on L reveals the existence of a characteristic length scale £at beyond which 
Ra(L) is enhanced up to unity despite the low interfacial transparency. The Andreev reflection 
enhancement is due to the strong mixing of particle and hole states that builds up in contacts ex- 
ceeding the coherence length £jv, leading to a small energy gap (minigap) in the density of states of 
the normal system. The role of the geometry of such hybrid contacts is discussed in the context of 
the experimental observation of zero-bias Andreev anomalies in the resistance of extended carbon 
nanotube/superconductor junctions in field effect transistor setups. 

PACS numbers: 74.45.+C, 74.50.+r, 73.23.Ad. 



I. INTRODUCTION 

The interest in proximity-induced supercon- 
ductivity in one-dimensional (ID) electron sys- 
tema 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ^ 10 ! 1 ^ 12 ! 13 has recently revived in 
the light of successful experiments on electron transport 
through carbon nanotubes contacted by superconduc- 
tors 1 ^* 1 ^* 1 ^ 7 ". Despite possibly strong electron-electron 
interactions 1 ^ 1 ^, in many situations transport prop- 
erties of metallic single-walled carbon nanotubes can 
be interpreted within a ballistic model assuming two 
conduction bands at the Fermi leveUS. Hence, supercon- 
ductor/carbon nanotube (S/CN) junctions can to some 
extent be viewed as an experimentally accessible case of 
ID ballistic proximity structures^ 5 -. 

Like in conventional normal metal/superconductor 
(N/S) junctions, the extent to which the proximity effect 
modifies the electronic properties of carbon nanotubes 
strongly depends on the quality of S/CN interfaces. In 
S/CN/S junctions with highly nontrivial end-bonding of 
the tubes it is possible to achieve high transparency con- 
tacts and observe induced supercurrents between the S 
banks^i 1 ^ 7 -. In a more conventional field-effect tran- 
sistor setup a superconductor is sputtered on top of a 
nanotube covering it from the ends and in this way con- 
necting it to the leadsi 5 -. Such contacts exhibit no ob- 
servable superconducting coupling, probably because of 
a Schottky barrier formed at the S/CN interfaces. Never- 
theless, in this case the proximity effect manifests itself as 
a pronounced zero-bias dip in the low-temperature resis- 
tance to which either of the S/CN interfaces contributes 
independently 15 . 

The sensitivity of the zero-bias resistance anomaly to 
the temperatur e) 15 ! 17 suggests that it can be attributed 
to the conversion of a normal current into a supercur- 
rent via the Andreev reflection process^ during which 
particles with energies much smaller than the supercon- 
ducting gap A are coherently scattered from an S/CN 



interface as Fermi sea holes back to the normal system. 
Under assumption of each of the S/CN interfaces act- 
ing independently^ 5 - and in the picture of non-interacting 
electrons, such an interpretation must reconcile with the 
well-established ID scattering model for a single N/S con- 
tact. However, for a point contact of average qual- 
ity (between metallic and tunnel regimes) the theory of 
Ref. |^ predicts a zero-bias resistance peak at temper- 
atures T < A/fe_Bj that is exactly the opposite to the 
experimental findings of Refs. ll5lll7l in the same temper- 
ature regime. 

Deviations of Andreev reflection physics in ID prox- 
imity structures from the standard model of Ref. l2ll 
have so far been ascribed to repulsive electron interac- 
tion^^LMJii or disorder in the normal channel^. In the 
present paper we show that clean non-interacting ID sys- 
tems can also exhibit unusual Andreev reflection proper- 
ties if the contact to the superconductor is not a point- 
like one. Such contacts naturally occur in field effect 
transistor setups due to a finite overlap between a nan- 
otube and a superconductor coupled in parallel. In par- 
ticular, in the device of Ref. Il5l this overlap was as large as 
lfj,m. To demonstrate the importance of the contact ge- 
ometry, we develop a scattering model for phase-coherent 
electron transport through a normal ID ballistic channel 
part of which is in parallel coupling to a 2D superconduc- 
tor via a low-transparency barrier (Fig. [Q. This model 
is in many aspects different from the device of Ref. 
and it is not expected to desribe all the experimental 
features. However, it captures the most essential, for 
our purposes, attribute of the S/CN contacts, namely 
their extended character. Moreover, the proposed geom- 
etry may serve as a minimum model accountin g for the 
zero-bias resistance features reported in Refs. I22l23l2 l] 
for extended planar contacts between ballistic 2D elec- 
tron systems and superconductors, whose cross-sectional 
structure is similar to that shown in Fig. ^ 

Our numerical simulation of elastic quasiparticle scat- 
tering shows that the probability of Andreev reflection 
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depends on the length L of the contact, approaching 
unity as L exceeds a certain length scale £jv larger than 
the coherence length in the superconductor £g. Most im- 
portantly, at zero energy high-probability Andreev reflec- 
tion occurs at any finite interfacial transparency for suffi- 
ciently long contacts. This is in sharp contrast to the sit- 
uation in point junctions^. To rationalize this result we 
perform a numerical analysis of the quasiparticle density 
of states (DOS) in the region of the ID system coupled 
to the superconductor in the limit L ^> £ N . The DOS is 
found to have a proximity-induced gap (minigap) at the 
Fermi level whose size E g k, (£,s/£n)A is much smaller 
than the gap A in the superconductor. The minigap 
Eg scales with the interfacial transparency, which implies 
that it is due to the formation of mixed particle-hole (An- 
dreev) states 2 ^. A comprehensive analysis of the energy 
dependence of electron scattering reveals that the gapped 
excitation spectrum in the proximity region results in the 
enhancement of the Andreev reflection probability Ra(^) 
at finite (but small) energies e < E g -C A followed by its 
decrease at intermediate energies E g < e < A. At the 
edge of the superconducting gap (e = A) the dependence 
Ra(s) exhibits one more peak typical for tunnel junc- 
tions^. These features dominate the bias voltage depen- 
dence of the differential resistance which at T ~ E g /ks 
displays a dip around the zero voltage similar to that 
observed in Ref. [Til 

Previously, zero-bias conductance anomalies 
have been extensively studied in mesoscopic 
superconducting contacts with semiconductors 

2fi.27.22.2S.29.3n.31.32.3.^34.23.24.3R.36.37 ^ ;||( , 

3^13,38,39,40,41,42,43,44 _ Thege gtudies have predominantly 
focused on the diffusive transport regime. According to 
the semiclassical scattering interpretation of Ref. I2H the 
excess conductance (i.e. exceeding the value predicted 
by the theory of Ref. I2H) is a signature of the correlated 
particle-hole motion arising from multiple Andreev 
reflections at the interface mediated by elastic scatterers 
in the normal system. Even for a low-transparency 
contact the cumulative Andreev reflection probability 
can be ~ 1 for trajectories hitting the interface many 
times provided that the area of the contact is sufficiently 
large. In the less explored regime of ballistic propaga- 
tion, a similar process, called sometimes reflectionless 
tunneling 29 , occurs in ballistic quantum wells in parallel 
long contacts with superconductors. In these setups, the 
multiple Andreev reflections are due to the back wall of 
the quantum well 2 V 5 i 35 . 

Therefore, the low-bias excess conductance discussed 
in the present paper is a pronounced case of the reflec- 
tionless tunneling in ballistic systems where Ra(c) can be 
interpreted as the cumulative Andreev reflection proba- 
bility due to the correlated quantum particle-hole motion 
in the region of the f D channel coupled to the supercon- 
ductor. It is also known that such correlations can lead to 
a minigap in the quasiparticle DOS 25 i 40 i 41 i 43 i 45 i 46 i 47 . We 
note that the previous studies of reflectionless transport 
Refs. I21] and |3f| dealt with multiple Andreev reflections 



semiclassically and at small energies |e| <C E g . Our quan- 
tum scattering approach is capable of describing the en- 
tire energy dependence of the subgap conductance which 
shows the crossover from reflectionless tunneling to in- 
dependent electron tunneling through the barrier at the 
N/S interface. Besides, our numerical technique allows 
us to tackle the realistic geometry of finite-length parallel 
N/S contacts and to obtain an accurate complete depen- 
dence Ra{L) which has not been studied in the previous 
models. 

The structure of the article is as follows. In Section ITU 
after a brief description of our system, we present the 
numerical results for the DOS in the ID channel. An an- 
alytical model is also developed that helps to rationalize 
the low-energy regime. The length and energy depen- 
dence of the Andreev reflection probability is analyzed 
in Section IlIII In Section llVI we summarize the implica- 
tions of our results with a concluding discussion on the 
bias voltage dependence of the resistance. 



II. TWO-GAP SPECTRAL PROPERTIES OF 
EXTENDED SUPERCONDUCTING TUNNEL 
CONTACTS 

In this section we study the density of states (DOS) in 
a quasi-one-dimensional electron system (Q1DES) cou- 
pled in parallel to a superconducting film via an inter- 
facial barrier. We consider the two geometries shown in 
Figs. Ola) and (b). The heterostructures are assumed 
two-dimensional and located in the plane x, z. Our re- 
sults can be easily extended to an out-of-plane periodic 
structure defining a quasi-two-dimensional electron sys- 
tem on the normal side. 



A. Description of the method 

To analyze the superconducting proximity effect in the 
QIDES we employ a numerical approach to solve the 
Bogolubov-de Gennes (BdG) equation 

(HA \fu(x,z)\_ (u(x,z)\ m 
\A -H* J \ v(x,z) )- e {v(x,z) ) 

for the electron u(x,z) and hole v(x,z) wavefunctions. 
The method allows for performing a straight-forward dis- 
cretisation on a real-space grid of the one-particle Hamil- 
tonian H = — (h 2 /2m)(d 2 + d^) + U(x 1 z)-fi, the pairing 
potential A = A(x, z), and the potential U(x,z) (to be 
defined later), fi and m are the chemical potential and 
the electron mass, both constant throughout the entire 
system. No translational invariance in the x-direction is 
invoked so that studies of the structures in Figs.^a) and 
(b) are possible. 
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FIG. 1: Quasi-one-dimensional electron system (Q1DES) 
coupled to: (a) semi-infinite and (b) finite length L supercon- 
ductor contacts. The interfacial barrier is depicted in black. 
The Q1DES width is chosen to be Wn ~ Af/2, where Af is 
the Fermi wavelength. 



The discrctized BdG equations read 

[e- (4 7 + C7(i,j) -n)]u(i,j) 
+7E;v "(*',/) ~ &(i,j)v(i,j) = 0, 



(2) 



[-e-(4r/ + U(i,j)-iM)]v(i,j) 

+ 7 E, y ^'j') + a(u>(m1 = o, 



where i and j refer to sites on a two-dimensional lattice in 
the x and z directions, respectively, and primes denote 
summation over nearest neighbours. The origin of the 
coordinate frame is indicated in Figs, da) andQJb) by 
the zero. If required, Eq. can be generalized to a 
position dependent effective mass for specific materials 48 . 

The parameters of the numerical scheme are as fol- 
lows. The potential U(i,j) is infinite everywhere outside 
the N and S systems. For every j within the materials, 
U(i,j) = for i < and U(i,j) = U Q > otherwise. A 
positive potential step U accounts for the fact that the 
coupling to the superconductor may result in a slight re- 
duction of the Fermi energy of the Q1DES in region (II) 
in Figs.ma) and^b) compared to that in the uncoupled 
region (I) (cf Ref. 1241) . This turns out to be important 
when considering scattering of quasiparticles incident at 
region (II), which is analyzed in the next section. The 
pairing potential in Eq. J2J) is assumed position- 

independent and equal to A • Sij (s-wave) in the super- 
conductor and zero everywhere else. Although the self- 
consistency is ignored, the stepwise order parameter has 
proved to be a satisfactory approximation for studying 
the proximity effect in clean systems 29-33 . The absolute 
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FIG. 2: DOS u(e) of a ballistic QlDES-superconductor sys- 
tem. The solid line corresponds to the semi-infinite geometry 
shown in Fig. 0^a) for which i/(e) is calculated directly from 
the Green function of the system. The dashed line corre- 
sponds to the case of a finite (but relatively long L/£s S> 1) 
S film [Fig. |T^b)] where we use the decimation technique. 



value of 7 is inverse proportional to the mesh param- 
eter a, which is varied until convergence of the results 
is reached. A sufficient condition is £s,Ai? 3> a, where 
£<? = hvp/2/S. is the superconducting coherence length 
and Xf (vf) is the Fermi wavelength (velocity). To sim- 
ulate the effect of a relatively thick superconducting film, 
we consider Ws/£,s — 15 for the spectral properties and 
Ws/£s — 50 for the results of Sect ion ITTTl Increasing this 
ratio does not have any quantitative effect at |e| < A. In 
particular, all features discussed below are already ob- 
served for Ws/£,s ~ 3 but with prominent finite-size ef- 
fects for high quasiparticle energies |e| > A. The width 
of the normal region is fixed to IFjv/Af ~ 1/2, allowing 
only one propagating mode along the Q1DES. The ratio 
£,s/Xf is chosen to be 2 (see also endnote 0) . 

A tunneling barrier at the N/S interface (dark area in 
Fig. ^) is introduced via an effective Hamiltonian equiv- 
alent to adding 



Et'-f' (tws - 7) feofy.i + 5j,i$f,o) u(i',f) = 0, 



E*v (7ivs - 7) {Sj,oSj',i + Sj,i8j<,o) v(i',f) = 



(3) 



to Eq. (|2J) for every i within region (II). Essentially, the 
above boundary conditions define the coupling between 
the normal and the superconducting systems via the in- 
terfacial constant jns- The latter may arise from a for- 
mal procedural that projects out the degrees of freedom 
within an insulating layer with a very high barrier, when 
neglecting the energy and momentum dependence of the 
penetration length. In what follows, we express all ener- 
gies in units of 7 for convenience. 
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FIG. 3: (a) Decomposition of the DOS u(e) to that aris- 
ing from the Q1DES (solid curve) and the superconductor 
(dashed curve). Both curves are scaled to compare with each 
other, (b) DOS in the proximity region for various interfacial 
transparencies. Inset: the minigap E g depends quadratically 
on the interfacial coupling jns at low transparencies, and 
starts deviating at higher values of jns- 



B. Quasiparticle density of states 

The proximity effect is reflected in the DOS, v(e), of 
the hybrid system plotted in Fig. for jns — 0.32 and 
U = 0. The solid line corresponds to the semi-infinite ge- 
ometry shown in Fig. ^a). In this case v(e) is calculated 
from the Green function G(i,j,i',j';e) of region (II) via 
the well-known relation v(e) — — ^ ^ X)j . G(i, j, i, j; e), 
where the summation is over all lattice sites in the hy- 
brid part of the junction; the factor N normalises the 
area under the curve to a reference unit and 3 means 
the imaginary part. The dashed line corresponds to the 
geometry of Fig.QJb) for a relatively long wire L/£s ^ 1- 
In this case, we obtain v(e) using a recursive technique 
(negative-factor counting 5 ?) that allows one to calculate 
the effective (renormalised) interaction between the nor- 
mal leads^S by projecting out the degrees of freedom of 
the middle region (II). 

Both approaches reveal the formation of a two-gap 
structure: a smaller gap (minigap) at E g rs 0.21A and 
the usual BCS singularity at e = A with a finite quasipar- 
ticle contribution at intermediate energies, E g < e < A. 
Figure Ola) shows the DOS separately for the Q1DES 
(solid line) and for the superconductor (dashed line) from 
which one can conclude that the smaller gap opens in 
the DOS of the Q1DES. This observation along with the 
dependence of E g on the coupling to the superconduc- 
tor shown in Fig. I^b) suggests that the minigap forma- 
tion is a signature of the superconducting correlations in- 
duced in the Q1DES. They are maintained in the course 
of multiple Andreev reflections in the channel which mix 
particle and hole states with energies below the effec- 
tive pairing energy coinciding with Eg&. At higher en- 
ergies Eg <C e < A the electrons and holes in the Q1DES 
are weakly correlated and hence can be treated as one- 
particle excitations. 

In mesoscopically large diffusive N/S systems the for- 
mation of the minigap has been stu died i n a number 
of theoretical papers (see, e.g, Refs. I4(l41ll4l . In the 
clean limit, the minigap structure has been analyzed to 
some extent in billiard geometries resembling quantum 
dota 45 i 46 . 47 . 

Below we develop an analytical model, close in spirit 
to our numerical approach, that provides a simple de- 
scription of the superconducting correlations in ballistic 
wires based on a ID BdG-like equation with an effective 
proximity-induced pairing energy E g . 



To calculate the DOS of the hybrid N/S system (re- 
gion (II) in Fig. ^) and, later, to study the scatter- 
ing and transport properties of quasiparticles incom- 
ing from normal region (I), we use a volatile numerical 
method used in studies of the magnetoresistance of hy- 
brid systems^iiS, phonon transport, and more recently 
in molecular electronic* 54 . It is based on recursive Green 
function techniques. Although some of the implementa- 
tion details may differ—, the main stages of the compu- 
tational scheme are explained in Ref. |52j. 



C. Proximity effect in a clean quantum wire: An 
analytical model 

Although the physical mechanism responsible for the 
minigap formation in ballistic 2D electron systems has 
been explored in Ref. |2^, the proposed method of deriva- 
tion of E g heavily relies on the following assumptions. 
Firstly, the pairing potential in the superconductor A(z) 
was assumed homogeneous. Secondly, a finite-thickness 
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normal system was modelled by a rectangular confining 
potential U c (z) enabling a plane wave description of the 
multiple reflection in the normal channel. The advantage 
of the model of Ref. |25| is that it allows one to obtain the 
minigap for an arbitrary interfacial transparency. In this 
subsection we present an alternative microscopic deriva- 
tion of Eg that does not rely on any particular models 
for A(z) and U c (z), but is restricted to low interfacial 
transparencies and low energies |e| ~ E g <C A. By fo- 
cusing on this case (weak-coupling regime) we would like 
to emphasize that the effects related to the minigap for- 
mation can be observed even in samples with average 
interfacial quality provided that the temperature is low 
enough. As in Ref. l25l we also assume the translational 
invariance along the N/S interface and different Fermi 
energies /is 3> (J-n and Fermi momenta ps S> pn on the 
S and N sides. 

It is convenient to rewrite the BdG equation for 
the two-component wavefunction ip p (z) = [u p (z),v p (z)] T 
in the superconductor (z > 0) in a more compact form: 



eo- 3 



2m 



d 2 z a - A(z)ia 2 



^(z)=0. (4) 



Here, p = p x is the momentum parallel to the interface; 
(72,3 and do are the Pauli and unity matrices, respectively. 
In the normal system (z < 0) the equation for ip p (z) is 



ecr 3 



+ ( E N + &dl - U c {zj) a Q ] M z ) = 0: (5) 



with Eff t s = fiN,s ~ p 2 /2m. The confining potential 
U c (z) defines a Q1D channel with a localized electron 
wavefunction cj)(z) in the z-direction. 

The interfacial barrier is assumed rectangular with the 
electron penetration length k,^ 1 = h/ftmU) 1 / 2 deter- 
mined by the barrier height U measured from the Fermi 
energy. Inside a high enough barrier one can neglect the 
energy and momentum dependence of the penetration 
length and write the BdG equation as [d 2 — Kfflj) p (z) — 
0, < z < a, where a is the barrier thickness. We intro- 
duce a special notation ipp( z ) f° r the BdG wavefunction 
inside the barrier to distinguish it from that outside the 
barrier. The continuity of the particle current imposes 
usual boundary conditions at the barrier walls, reading 

VV(0) = W(0), Ma) =Vp(o), (6) 
9 Z -0 P (O) = <9 2 Vp(0), <9j^ p (a) = d z tp p (a). ( 7 ) 

The solution inside the barrier satisfying the boundary 
condition® is = Vfe(0) + ffi%g 

Inserting it into the boundary conditions Q for the 
derivatives, we have: 



0,1^(0) 
d z if) p ( a ) 



Kl/l p ( a ) 



K t ip p (a), 

-K t 1pp(0), 



(8) 
(9) 



where n = Kocotanhft a and n t = kq/ sinhK a. Equa- 
tions iJSJ and © serve now as effective boundary condi- 
tions for the BdG equations in the superconductor and 



the normal system. In the limit sinhKoa — > oo, the cou- 
pling between the "normal" and the "superconducting" 
functions vanishes, which is described by Eqs. (jHJ and 
@ with zero right-hand sides. 

We use boundary conditions J8J and © to describe 
Andreev reflection at the superconductor-QIDES inter- 
face under the assumption that the influence of the 
Q1DES on the superconductor can be neglected. To pro- 
ceed, it is convenient to include the boundary condition 
@ into the BdG equation 10} by introducing appropriate 
delta-function terms as follows 



6(73 +[E S + ^dj + Us(z) cto - A(z)ia 2 xp p {z) 



- K t h S(z 

2m "V" 



a)Vv(O). 



(10) 



We note that the admitted singular potential Us{z) = 

■§^S(z — a)(d z — k) reproduces Eq.|5|with zero right-hand 
side ("isolated superconductor"). 

The penetration of Andreev bound states into the su- 
perconductor at low energies is described by a particular 
solution of Eq. lfTU|) generated by the right-hand side con- 
taining the "normal" function tp p (0). It can be expressed 
in terms of the matrix Green function of Eq. (|10[1 whose 
matrix elements are constructed from the quasiparticle 
G Pje (z,z') and condensate (Gorkov's) F Pj€ (z, z') Green 
functions, namely, 



ip p (z) 



= _Kth 2 _ ( G p , e (z,a) -F_ p _ e (z,a) 
2m \F PjC (z,a) G„ Pj _ e (z,a) 



V>p(0).(ll) 



Here the Green functions satisfy boundary condition @ 
with zero right-hand side. Inserting this solution into the 
boundary condition © at the " normal" side and neglect- 
ing both energy and momentum dependence of the Green 
functions under conditions |e| -C A and p < p^ *C ps, 
one finds 



0^p(O) + k^p(0) = -l—Fw 2 M°)> 
2m 



(12) 



where F = F p — e —o{a, a) is the condensate Green function 
taken at the boundary of the superconductor. We have 
omitted the terms proportional to G since for |e| <C A 
they would only result in a shift of the dispersion. 

The right-hand side of the boundary condition Ql'ifl. 
which is off-diagonal in the particle-hole space, takes into 
account the conversion of a particle into a hole (and vice 
versa) due to Andreev reflection, that occurs simultane- 
ously with normal scattering. In a narrow quantum wire, 
whose thickness is of order of the Fermi wavelength, the 
anomalous term in the boundary condition l|12|) gives rise 
to an effective pairing energy between particles and holes 
in the wire. Indeed, combining the equation of motion 
(JSJ and the boundary condition (|12f> . one can write 



eer 3 



E 



n-i)- 



N 



U c (z) + U N (z] 



CO 



(13) 



where the singular potential Un(z) = — ^<5(z)(9 z + k) 
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is equivalent to the boundary condition 1)12(1 with zero 
right-hand side. For a weakly coupled Q1DES, the spa- 
tial dependence of the BdG function ipp(z) w tp p (j)(z) is 
almost unaffected by tunneling. Therefore, multiplying 
Eq. ED by 4>(z) and integrating over z, one obtains the 
following one-dimensional equation: 



ecr 3 



( p%-p 2 \ 

\ 2m J 



(T - E q ia 2 



E g EE (k 4 /i 2 0(O)/2to) F. 



= 0, (14) 
(15) 



Egioi plays the role of the effective singlet pairing en- 
ergy in the wire; pf denotes the Fermi momentum in the 
Q1DES. 

According to Eq. (|14fl . the excitation spectrum in 
the Q1DES is ef = ±[v F (\p\ - pf) 2 + Eg} 1 / 2 with the 
Fermi velocity vp = pFjm. It has an energy gap given 
by Eq. 1(151) and, hence, the DOS of the normal sys- 
tem displays a BCS-like singularity at E g . To esti- 
mate Eg, one can use the condensate Green function 
of a superconductor with a homogeneous pairing po- 
tential A at zero energy and parallel momentum, F « 
WV'Ep, A/(A 2 +v 2 s ( Pz - Ps ) 2 ), where v s = Ps /m. The 
integration over p z gives F ks 1/hvs- The boundary value 
(f)(0) of the transverse function can be estimated using the 
unperturbed boundary condition 0(0) = — K~ 1 d z <fi(0), 
where on the right-hand side one can use the " hard wall" 
wavefunction 4>(z) = (2/Wn) 1 ^ 2 svhttz/Wn, which gives 
|0(O)| w h- 1 {2/W n ) 1 ' 2 {-k/W n ). Thus, the effective pair- 
ing energy is 
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FIG. 4: Andreev reflection probability for various N/S cou- 
pling constants as a function of the length of the proximity 
region (II) in Fig. |5fb). Ra scales with the interfacial trans- 
parency and at short lengths is quadratic in L. For L much 
larger than the proximity-induced coherence length £jv fa 5£s, 
the probability Ra reaches unity as a manifestation of the re- 
flectionless tunneling despite the low interfacial transparency 
(see inset for fMS = 0.2). 



III. QUASIPARTICLE SCATTERING: LENGTH 
AND ENERGY DEPENDENCE 



E B=rxT—-rA E o, (16) 

WnPS smh n a 

with Eq = h 2 ir 2 /2mW N being the energy of the lowest 
occupied subband in the quantum well. Equation 1161) 
is equivalent to the one obtained in Ref. |2^ for a strong 
delta-shaped barrier. 

Equation (|15f) for the minigap E g provides a link to 
the numerical approach and results discussed earlier. Ac- 
cording to Eq. I |l 5(1 . the size of the minigap depends on 
the parameter K t ti 2 /2m that characterizes "hopping" be- 
tween the systems (see Eq. (111(1 ). This parameter rep- 
resents a direct analogue of the coupling constant "/ns 
that determines the size of the minigap in the DOS in 
our numerical study [see Fig. [Sfb)] . Since E g is quadratic 
in Kfh 2 /2m, the numerical value of the minigap should 
scale with jns a s 

E g oc Ins, (17) 

which can indeed be verified numerically (sec inset in 
Fig. GJb)). In the next section we will see that the 
parabolic dependence of the effective pairing energy l(17|) 
on 7jv5 can also be extracted from calculations of the 
Andreev scattering probability. 



In this section we discuss electron scattering properties 
that can be used as an independent and more complete 
probe of the proximity effect in finite length parallel N/S 
contacts. For definitcness we consider particles in the left 
region (I) of Fig. ^ propagating to the right. When in- 
cident at the boundary with the proximity region, these 
may be: (a) Andreev reflected, namely, converted into 
outgoing holes with the probability Ra, (b) normally re- 
flected as outgoing particles, i.e., without Andreev con- 
version, with the probability i?jv, and (c) normally trans- 
mitted as particles with the probability Tjv either in the 
region (II) of Fig. |Ja) or in the right region (I) of Fig. 
^b). Finally, the probability of being Andreev transmit- 
ted to the right is determined via particle conservation, 
namely, 1 — Ra — Rn — Tjy. In our calculations, this is 
ensured by the unitarity of the scattering matrix. 

We examine first the dependence of the zero-energy 
Andreev reflection coefficient Ra{^ — 0) on the length 
L of the proximity region (II) in Fig. ^b) for different 
values of the coupling parameter jns and without any 
potential mismatch at the (I)/ (II) boundary (U = 0). In 
conventional N/S/N structures the Andreev coefficient is 
known to scale as Ra ~ (A L) 2 for L much shorter than 
the coherence length £,g^- According to the results of the 
previous section, in our case the effective pairing energy 



Eg [Eq. djl] should act as A and therefore we expect 
that Ra ~ {EgL) 2 or, according to Eq. El Ra ~ Ins^ 2 
for short enough contacts. This scaling is demonstrated 
in Fig. 01 by the convergence of the appropriately nor- 
malized Ra curves and their parabolic shape at short 
lengths. As shown in the inset, there is a characteristic 
length £jv ~ 5£s beyond which the Andreev probability 
Ra{L) approaches its unit limit. Moreover, the ratio of 
£,n/£s coincides with the ratio of the gaps A/E g sa 5 
found from the analysis of the DOS in the previous sec- 
tion: 



6V/£ S = A/Eg 



(18) 



The overall length dependence implies that the reflection- 
less tunneling builds up due to the strong mixing of par- 
ticles and holes in long channels. In p articular, the semi- 
classical approaches of Ref . l27l andl35l interpret reflection- 
less tunneling in terms of the increase in the cumulative 
Andreev reflection probability with increasing number of 
single Andreev reflections at the N/S boundary in the 
limiting case of an infinitely long interface £/£_/v - * oo. 

We now turn to the discussion of scattering of finite- 
energy quasiparticles in the semi-infinite geometry of the 
proximity region [Fig. [Ja)] where the reflectionless tun- 
neling is most pronounced. We also take into account 
a finite potential step U at the boundary between the 
normal (I) and proximity (II) regions that, as has been al- 
ready mentioned, may arise from the modification of the 
Fermi energy in region (II) due to the coupling to the 
superconductor. The energy dependence of the Andreev 
reflection coefficient is plotted in Fig. OJa) for various U 
and 7jvs = 0.32. At low energies e < E g = 0.21A the 
shape of the dependence Ra(c) resembles that of high- 
transparency N/S point contacts discussed by Blonder, 
Tinkham and Klapwijk (BTK)2i. If there is no poten- 
tial step U between the normal (I) and proximity (II) 
regions, the probability Ra equals unity and starts to 
drop at e > E g . For U a ^ 0, finite normal reflection 
builds up (shown in Fig. EJa)) which results in smaller 
zero-energy values of Ra- The appearance of the second 
narrow peak at e — A 58 manifests the crossover from 
the reflectionless tunneling regime, which involves a two- 
particle process, to the usual independent electron tun- 
neling through the barrier. 

In Fig. Efb) we demonstrate that the low- and high- 
energy peaks in the dependence of Ra can be indepen- 
dently fitted by the BTK model^I. To fit the low-energy 
behaviour we use the formulas of Table II in Ref. with 
A btk = Eg and a small barrier parameter Z — 0.278. 
For the tunneling peak we use the same formulas with 
Abtk = A and the large barrier parameter Z = 28. In 
either of the above limiting cases, the fit is almost per- 
fect. To describe the crossover between them, a more 
general analytical model is needed. 

In Fig. all non-vanishing scattering coefficients are 
plotted for: (a) U — 0.8 and (b) no barrier between 
the normal (I) and proximity (II) regions. For U =/= 
there is normal reflection of particles at the (I) / (II) in- 
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FIG. 5: (a) Andreev reflection coefficient for the geometry of 
Fig. H^a) with ^ns = 0.32 and various values of the potential 
Uo at (I) / (II) boundary. For U — the Andreev probability 
is exactly 1 for energies below Eg. (b) To fit the low- and 
high-energy peaks we use the formulas of the BTK model— 
with parameters Abtk = E g , Z = 0.278 and Abtk = A, 
Z = 28, respectively. 



terface caused by the potential mismatch. For U = it 
vanishes not only below the minigap e < E g but also at 
the intermediate energies E g < e < A. This is due to 
the specific geometry of our tunnel junction where quasi- 
particles with intermediate energies are mainly transmit- 
ted through the channel experiencing low-probability An- 
dreev reflection (cf the behaviour of T/v(e) and i?^(e)). 
At higher energies e > A, when the superconductor 
becomes transparent for quasiparticles, the lack of the 
translational invariance of our system causes considerable 
normal scattering and oscillations of all the coefficients 
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FIG. 6: Andreev reflection Ra, normal transmission Tn, and 
normal reflection Rn coefficients for the geometry of Fig. 0^a) 
with 7jvs = 0.32 and various potential landscapes in region 
(II). 



due to the finite thickness of the superconductor. 



IV. CONCLUDING REMARKS 

Scattering coefficients are related to the two-, three-, 
and in general multi-probe conductances of N/S sys- 
tems (see, e.g., Refs. l2lL l33f) . Therefore, the proximity 
effect discussed in the previous section should be observ- 
able in measurements of the current-voltage characteris- 
tics of such hybrids. In order to explore this possibility 
and to elaborate on our discussion of the experimental 
reportsi£i22i22i24, we conclude by focusing on the two- 
probe differential conductance g(eV) = dl/dV of the 




eV/A 



FIG. 7: Differential conductance of the N/S hybrid in 
Fig. IE[&) at ultralow temperature T <C Eg/ks- 



semi- infinite geometry of Fig. HJa) . This is given bj^ 



9(eV) = 




where fP^ = {cxp[(e =p eV)/k B T] + l}" 1 with (-) for 
(p)articles and (+) for (h)oles. The bias energy eV is 
introduced as the difference between the chemical poten- 
tials in the normal region and in the hybrid part of the 
junction, with the latter taken as reference. 

At zero-temperature, Eq. I|19ll reduces to g(eV) = 
1 + R p A (eV) - R p N (eV) = 1 + R A (-eV) - R%(-eV). 
Hence, for a small barrier between the normal (I) and 
proximity (II) regions of the wire, the dependence of 
g(eV) at ultralow T reflects mainly the energy depen- 
dence of the Andreev probability [Fig. EJa)]. The same 
is true for the length dependence of the zero-bias con- 
ductance. In Fig. we plot the differential conductance 
for several values of interfacial coupling 7^,5 (i.e., Eg) at 
T = O.OlA/fcs. Unlike the tunneling peaks at ±A/e, 
the proximity-induced anomalies at the minigap energy 
±E g /e exhibit a strong dependence on jns (cf Fig.[3fb)). 

In Fig. IS] the evolution of the differential resistance, 
which is defined as the inverse of Eq. ^3 is shown as a 
function of temperature. At intermediate T ~ Eg/ks, 
features at the scale of the minigap are smeared and 
the resistance exhibits an overall dip as a result of the 
reflectionless tunnelling. With decreasing the temper- 
ature to T < Eg/ks the resistance curve develops a 
finer structure reflecting the energy dependence of the 
Andreev reflection probability. For vanishing potential 
step Uq (Fig. |SJa)), there are two minima symmetric to 
zero bias at the energies of the superconducting gap. In 
addition to those, for Uq 5^ (Fig. 02b)) the finite normal 
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FIG. 8: Evolution of the differential resistance of the N/S 
hybrid in Fig.^a), as a function of temperature from ksT ~ 
Eg < A to T <C -E 9 for (a) vanishing and (b) finite potential 
step between regions (I) and (II). 



reflection at |e| < E g leads to a zero-bias resistance peak 
superimposed on the Andreev dip. 

According to Fig. 0] for pronounced Andreev reflection 
the contact length L must be of order of the proximity- 
induced coherence length £jv = 5£s- For O.Ol^m < £g < 
0.04/im typical for Nb electrods, £jv can be estimated as 
0.05/im < £jv < 0.2//m. On the other hand, in the device 
of Rcf. |l5| the contact overlap between the nanotube and 
the superconductor was rather large, about 1/im. There- 
fore, the condition L > £/v could be met leading to the 
observed zero-bias reduction of the contact resistance. 
Another feature of our dV/dI(V) curves, namely the ap- 
pearance of a small zero-bias peak superimposed on the 
Andreev dip due to finite normal scattering at very low 



T (Fig. IHfb)), is also consistent with the experimental 
findings. 

We note that in the experiment of Ref. the low-bias 
behaviour of the resistance was sensitive to a gate volt- 
age applied to the carbon nanotube. In our model the 
effect of the gate voltage can be incorporated into the 
difference U a between the Fermi energies in the normal 
(I) and proximity (II) regions in Fig. 2] We have focused 
on the most interesting case of relatively small U when 
normal scattering does not impede the conversion of the 
quasiparticle current into the supercurrent in the prox- 
imity region. As U increases, the low-bias resistance 
dip in Fig. |5Tb) eventually evolves into an overall peak 
above the normal state value like in non-ideal N/S point 
contacts2i*S. 

We emphasize that the zero-bias anomaly discussed 
here is a property of a single parallel N/S contact. Al- 
though as argued in Ref. Ha in their S/CN/S devices the 
two CN/S interfaces acted independently, the role of the 
interelectrode coupling remains unclear. Such a ques- 
tion has been investigated numerically in Ref. for a 
somewhat simpler system where a carbon nanotube is 
connected to a normal metal and a superconductor via 
tunnel barriers (N/CN/S). It was shown that resonant 
tunneling through Andreev levels in the nanotube can 
significantly increase the low-bias subgap conductance 
similar to the situation in mesoscopic N/quantum dot/S 
structures^. For a comprehensive theory of transport 
in S/CN/S hybrids this aspect together with the contact 
geometry and electron interaction effects must be taken 
into account. In addition, one should bear in mind that 
in the experimental realizations a number of CNs have 
been contacted in parallel. 

In Refs. and |3 a strong zero-bias suppression of 
the resistance was found in ballistic 2D electron sys- 
tems in extended planar coupling to superconductors 
at T < A/ks- These systems can be considered as a 
generalization of that shown in Fig. ^ For perfect pla- 
nar interfaces, individual channels with possibly different 
barriers and interfacial transparencies 59 add up indepen- 
dently. However, interchanncl mixing must be considered 
for rough surfaces. The same applies when considering 
experiments in quantum wires^* with few propagating 
modes. It is worth noting that in this case the behaviour 
similar to the low-T differential conductance of Fig. [7| 
was observed. We believe that the transport anomalies 
observed in lDi£, quasi-lD^ and 2D22i24 systems have 
a contribution of a common nature stemming from the 
proximity-induced mixing of particles and holes which 
mediate the conversion of a normal current into a super- 
current along the contact on the scale of the coherence 
length £tv = £,sA/E g and at energies smaller than the 
minigap E g < A. 
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